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ABSTRACT 

Context. Massive galaxy clusters at intermediate redshifts act as gravitational lenses that can magnify supernovae (SNe) occurring in 
background galaxies. 

Aims. We assess the possibility to use lensed SNe to put constraints on the mass models of galaxy clusters and the Hubble parameter 
at high redshift. 

Methods. Due to the standard candle nature of Type la supernovae (SNe la), observational information on the lensing magnification 
from an intervening galaxy cluster can be used to constrain the model for the cluster mass distribution. A statistical analysis using 
parametric cluster models was performed to investigate the possible improvements from lensed SNe la for the accurately modeled 
galaxy cluster A1689 and the less well constrained cluster A2204. Time delay measurements obtained from SNe lensed by accurately 
modeled galaxy clusters can be used to measure the Hubble parameter. For a survey of A 1689 we estimate the expected rate of de- 
tectable SNe la and of multiply imaged SNe. 

Results. The velocity dispersion and core radius of the main cluster potential show strong correlations with the predicted magnifica- 
tions and can therefore be constrained by observations of SNe la in background galaxies. This technique proves especially powerful 
for galaxy clusters with only few known multiple image systems. The main uncertainty for measurements of the Hubble parameter 
from the time delay of strongly lensed SNe is due to cluster model uncertainties. For the extremely well modeled cluster A1689, a 
single time delay measurement could be used to determine the Hubble parameter with a precision of ~ 10%. 

Conclusions. Observations of SNe la behind galaxy clusters can be used to improve the mass modeling of the large scale compo- 
nent of galaxy clusters and thus the distribution of dark matter. Time delays from SNe strongly lensed by accurately modeled galaxy 
clusters can be used to measure the Hubble constant at high redshifts. 

Key words, cosmology: gravitational lensing - supernovae: general - galaxies: clusters: general - galaxies: halos - dark matter ~ 
cosmological parameters 

1 . Introduction genitor systems of SNe and, if the cluster potential can be esti- 

. mated properly, extend the SN la Hubble diagram up to redshift 

^ Massive clusters have been successfully employed as grav- ^ ^ 3 possibly even higher. Such measurements would dra- 
_Cd , itational telescopes (also known as Zwicky telescopes) m ^atically extend the redshift baseline for which the expansion 
order to probe astronomical objects at very high redshifts history of the universe can be probed, which could turn out to be 
(e.g.,'Kneib et al.ll2004|; [ Bradley et al. | |2008| ; | Bradac et al. || 2009t ^^^^^^-^^ for understanding the nature of dai'k energy. 
[Zheng et al. 2009). More recently, the use of near-IR observa 
tions to look for gravitationally magnified supernovae (SNe) be 
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A first transient object behind the massive cluster A1689 was 



... . f , , f • • J n . J found in a ground based pilot survey at ESO using the ISAAC 

hind massive galaxy clusters has been investigated. Such data fr: — V , . " ir>r.ok ^nrr — ■ . — ^ 1 ll -irvrvnL 

. . , . ° , ^ , • , , , ..c . (Moorwood et al. 1998) camera on VLT ( Stanishev et alJI2009t 

give insight into the star formation rate at high redshifts, the pro- . ^ , onnn u n t j n tt V- — 

° ° to ' r Goobar et al. 2009, hereafter Paper I and Paper II, respectively). 

~ , ^ . ITtt t The transient was consistent with a reddened Type IIP SN at z 

Send ojfpnnt requests to: T. Riehm _ q y^^j^jj ^ lensin ma nification Am - 1 4 ma Other SN 

* Based on observations made with the NASA/ESA Hubble Space ,■ , , ^ ensing magni ca ion m . mag. 

Telescope, obtained from the Data Archive at the Space Telescope candidates have been found in a surve>; with the newer HAWK-I 

Science Instimte, which is operated by the Association of Universities (Pirard et alj|2004; Casali et al. 2006; Kissler-Patig et al.||2008D 

for Research in Astronomy, Inc., under NASA contract NAS 5-26555. camera at VLT ( Amanullah et al, 201 1). Although not optimally 

These observations are associated with programs # 9134, 9289 and cadenced for the purpos e, HST detections o f lensed SNe may 

10150. also result from CLASH dPostman et al.l201lh . part of the Multi- 
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Cycle Treasury (MCT) program targeting 25 massive cluster, 
both at optical and near-IR wavelengths. 

In this paper we investigate the impact observations of gravi- 
tationally lensed SNe could have on the modeling of galaxy clus- 
ter potentials. We will focus on the detection feasibility with ei- 
ther an 8-m class NIR survey like the one with HAWK-I (FOV 
» 56 arcmin^), or with the Wide Field Camera 3 instrument 
on board the Hubble Space Telescope (HST/WFC3) with a 10 
times smaller FOV. Since SNe la have small intrinsic luminosity 
dispersion (after correction for lightcurve shape and reddening), 
we can estimate the absolute magnification of the SNe and thus 
break the so called mass sheet degeneracy of gravitational lenses. 
This degeneracy implies that the density distribution of the lens 
can always be rescaled and a constant-density mass-sheet added 
such that the, also properly rescaled, source plane is projected 
onto the same observed images. Thereby, by breaking this de- 
generacy, direct constraints on the distribution of the dark matter 
component in galaxy clusters can be obtained. 

If a SN is strongly lensed, a measurement of the time de- 
lay between the transient in the multiple images could poten- 
tially constrain the Hubble parameter (Refsdal 196 4j), and thus 
dark matter and d ark ener gy parameters (|Goobar et al.l 120021 
iMortsell & Suness on 2006; ISuvu et all l2010h . Though not as 
precise as other cosmological tests to study dark energy, the time 
delay technique has the major advantage of measuring cosmo- 
logical parameters at redshifts where few other probes are cur- 
rently available. Given the transient nature of SNe, the time de- 
lay between multiple images could potentially be measured to 
very high precision. The main limitation of this technique is the 
strong degeneracy between the lens mass mode l and t he Hubble 
parameter. Ht\ (e.g..[Wa mbsganss & Paczvnski 1 19941 : IWitt et alJ 
120001: lKochanekll2003 iZhao & Oini,2003.) . However, observing 
cluster lenses with potentials well constrained by a large number 
of already known multiple images or, in the case of a strongly 
lensed SNe la, direct mag nification informatio n, this degener- 
acy can be diminished (Q guri & Kawandl2003h . This technique 
is applicable to any object with a variable light curve which 
would make it possible to determine a time delay between dif- 
ferent images, including quasars (QSOs) and gamma-ray bursts 
(GRBs). Although, there are ~ 100 strongly lensed QSOs cur- 
rently known, time delays have been measured for only ~ 20 of 
them ('O gurill2b07l) . This is due to the fact that most QSOs show 
variability only at a fairly low level, making time delay deter- 
minations challenging. Time delay measurements of QSOs be- 
hind cluster lenses have so far onl y been poss ible in two cases: 
SDSSJ1029+2623 (three images. Ilnada et all [2 006) and SDSS 
J1004-H41 12 (five images. iFohlm eister et al.' 2008). Image mag- 
nifications of about Am ~ 4.5 mag have been recorded. However, 
the mass distributions in these clusters are not modeled well 
enough to constrain the Hubble constant. Although the time de- 
lay of a GRB, could be determined with very high precision, they 
are expected to be too rare to admit a survey of lensed GRBs 
behind clusters. Therefore, in this paper we will focus on time 
delays for multiply imaged SNe. 

The paper is structured as follows: In Sect.|2] we describe the 
cluster model used in our analysis. Sect.[3]investigates the obser- 
vational prospects for lensed SNe behind a cluster Constraints 
on the cluster modeling and the Hubble constant from lensed 
SNe are studied in Sects. |4]and|5] Finally, we discuss our results 
and conclude in Sect. |6] 

Throughout the paper, we assume the cosmological parame- 
ters to be Q.M = 0.3, Q.A = 0.7, and //q = 70 km s"' Mpc"'. 



Fig. 1: Magnification map of A1689 overlaid on top of an HST 
ACS image of the cluster. The Am = 1,2,3 mag (pink, blue, 
magenta) contours are shown for a source at z = 2. The green 
(orange) circles indicate the positions of strongly lensed galax- 
ies with spectroscopic redshifts and time delays below (above) 
5 years. Red points indicate strongly lensed sources with photo- 
metric redshifts. 




197.85 



m.875 
R.A. (2000.0) 



2. Cluster modeling 

The focus of this paper is the massive cluster A 1689 at redshift 
z = 0.187 which is one of the best studied clusters. For this clus- 
ter there are 34 known multiply imaged background galaxies and 
a total of 114 images. For 24 of the 34 background systems, there 
are secure spectroscopic redshifts available, ranging fr om z = 
1.1 to 4.9 (iBroadhurst et al.ll2005l: [Limousin et alj|2007h includ- 
ing a total of 82 images (see Fig.[Tl). For this cluster, a detailed 
mass model is available consisting of 272 parametrized pseudo- 
isothermal elliptical mass distributions (PIEMD ). This is an up- 
dated version of the mass model presented in iLimousin et al.l 
(l2007h . Each potential is characterized using seven parameters: 
the center position (RA, DEC), the position angle 6, and the pa- 
rameters of the mass profile: the central velocity dispersion cr, 
the ellipticity e, the core radius rcoie and the cut radius r^m- The 
model has 33 free parameters describing the mass distribution 
of the two large scale dark matter clumps (Clump 1 and Clump 
2), the dark matter halos of three individual galaxies which were 
found to play an essential role for producing multiply imaged 
systems (BCG, Galaxy 1 and Galaxy 2), as well as two scal- 
ing relations based on luminosity for the remaining identified 
cluster galaxies, L'r^ut and L*cr (compare Table [TJ. The latter 
parameters are used to assign masses to the cluster galaxies, 
which are mainly early-type galaxies, via the scaling relations, 
rcut(L) = reut(L/L*)'^^ and cr(L) = a-(LIUyi\ where U is the 
luminosity of a typical galaxy in the cluster 

For investigating the power of lensed Type la SNe observa- 
tions for constraining cluster models we are also studying cluster 
A2204 at redshift z - 0.1524. For this cluster only one multi- 
ply lensed image system at z = 1 .06 is known, making it much 
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Table 1: Parameter values inferred for the dark matter clumps considered in the optimization procedure for galaxy cluster A1689. 
Coordinates are given in arcseconds with respect to the BGC. Error bars coiTespond to Icr confidence level as infeiTed from the 
MCMC optimization. The model has a total of 33 free parameters. Values in brackets are not optimized. 

Clump RA DEC e 6 ''core (kpc) r^^ (kpc) cr (km s"') 

Clump 1 Ol^^ -8.5!;;:* OlF™ 90.4!|;f WLT^q [1515.7] lU5.5t_{1j 



Clump 2 -70.7!|^ 49.r;5° 0.77!!!!!* 78.6!;* 67.9!^^ [501.0] 647.rf^ 

BCG -1.1^!!^ O.r?* 0.48+™ 66.0i| 5.3+M ISO.S!"^ 441.7!|i-^ 

■1 + 13:3 OAT Q + tif 



'-0.2 -^-0.4 --"-m --"-_6,3^ ^^-V^y^ 



Galaxy 1 [49.1] [31.5] 0.70!»- » 113.6!;^;^ 25.6!^;; i(>i-2 2ii 262.8 



Galaxy 2 -^5J*_°-*^ 31.6':ll 0.80^™ 45.5!;''f 17.6!|-^ 182.8!|-f 425.6 

L' elliptical galaxy ... ... ..." . . . "" [0.15] 54.5!^^ 154.6!^^ 

Table 2: Parameter values inferred for the dark matter clump and the scaling relation considered in the optimization procedure for 
galaxy cluster A2204. Coordinates are given in arcseconds ■with respect to the BGC. Error bars correspond to Icr confidence level 
as inferred from the MCMC optimization. The model has a total of 5 free parameters. Values in brackets are not optimized. 
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less constrained than A1689. The mass model for this cluster 
consists of one main potential, describing the smooth cluster po- 
tential, and 34 cluster member galaxy potentials. All potentials 
are described by PIEMD profiles and the model has a total of 5 
free parameters: the ellipticity, e, position angle, 6, core radius, 
'"core, and velocity dispersion, cr, of the main potential as ■well as 
the mass-luminosity scaling relat ion for the cluster ga laxies, L*cr 
(compare Table |2] for details see iRichard et al Mm . 

LENSTOOL is a soft^ware package for modeling 
mass distributions of galaxies and clusters in the 
strong and ■weak lensing regime (iJuUo et all 120071 
http://www.oamp.fr/cosmology/lenstool/). Monte Carlo Markov 
Chains (MCMC) are used to constrain the free parameters of the 
cluster model derived from observational data of the background 
galaxies. The output contains a chain sampling the probability 
distribution of parameter values, each of which corresponds to 
a specific model realization. The density of parameter values 
gives a measure of the corresponding probability distribution 
for the parameters and can be used to estimate the errors of 
the corresponding parameters. For each of the realizations, the 
magnification and time delay function at any given position 
behind the cluster, can be computed. 

Assuming that we observe the magnification or time delay 
of a multiply imaged SN, we are able to rule out all realizations 
that do not agree with these (simulated) additional constraints. 
By comparing the remaining realizations to the original chain, 
we are then able to judge how powerful the additional constraints 
are in constraining the cluster potential. 

3. Supernovae behind clusters 

To investigate how efficient massive galaxy clusters are for de- 
tecting distant SNe and putting constraints on cluster modeling 
or cosmology, a good estimate on how many SNe behind the 
cluster that can be expected during a given survey time is needed. 
For example, for a 5 year monthly survey of one very massive 
A1689-Uke cluster with HAWK-I, it has been shown that the to- 
tal number of SNe expected in the background galaxies is on 
the order of 40-70 SNe, out of which approximately 20-30 are 
SNe la (depending on the underlying rates estimate for the vari- 
ous SN types, see Paper II). 

For the purpose of studying time delays using galaxy cluster 
A 1689, we focus primarily on potential SN explosions in the 24 



known multiply lensed background galaxies with spectroscopic 
redshifts. Ten multiply imaged galaxies with photometric red- 
shifts only are not taken into consideration to avoid additional 
sources of error due to uncertainties in photometric redshifts. 
Furthermore, with more data being collected, new multiple im- 
age systems might be detected. The numbers we present below 
should therefore be considered as lower limits. 

Since the rate of SNe is expected to be coupled to the star- 
formation rate (SFR) in these galaxies, we use rate predictions 
for the strongly lensed systems derived from local estimates of 
the SFR. As the UV luminosity is dominated by the most short- 
lived stars, it is closely related to star formation. Thus, we use 
Lasoo, the flux at rest-frame A^ff - 2800 A (redshifted to opti- 
cal bands) as a tracer of the SFR in the strongly lensed galaxies. 
Since there are very deep HST observations of A 1689 with mul- 
tiple filters, these estimates are remarkably precise. 

The absolute L2800 magnitude is derived after taking into ac- 
count the distance modulus, K-corrections, extinction and lens- 
ing magnification from the cluster model. Finally, we use the re- 
lation between L2800 and SFR from .Dahlen et al. (2007) to relate 
the flux to star formation, 

SFR(Moyr-') = L2800 (erg ■ s-'Hz-') x (7.0 • lO^^)"' . (1) 

The expected rate for core collapse SNe is calculated from 

'•cc(z) = kl^ X SFR(z) , (2) 

where kf = 0.007Mq' is estimated using a Salpeter IMF and 
a progenitor mass range of between 8 and 50 solar masses. For 
SNe la, we assume that the rate has two components, one propor- 
tional to the S FR and one proportional to stella r mass, according 
to the model of lScannapieco & Bildstenl(l2005l) . The rate is given 
by 

ria(f) = A ■ SFR(0 + B ■ M{t) . (3) 

The mass, M, of i ndividual galaxie s are calculated using mass- 
to-light ratios from lBell et al.l (l2003h . The resulting rates for core 
collapse SNe and SNe la in multiply lensed background galaxies 
with spectroscopic redshifts and predicted time delays below 5 
years are summarized in Table [3] In total, we expect ~ 0.6 core 
collapse SNe and ~ 0.05 SNe la per year to explode in the known 
multiply imaged galaxies with spectroscopic redshifts. 
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Table 3: All image pairs of multiply lensed background galaxies with spectroscopic redshift and a predicted time delay of less than 
5 years. Column 3 gives the predicted time delay between the images. A negative time-delay indicates that image 2 will occur 
before image 1. Columns 4 and 5 give the predicted magnification for image 1 and image 2, respectively. In column 6 and 7 our 
rate estimates are shown for core collapse SNe and SNe la, respectively. Columns 8-14 indicate the detectability of a lensed SN 
for the different SN types assuming a limiting magnitude of 25.25 mag (Vega) in F140W, e.g. in a survey using 1/3 orbit with the 
HST/WFC3 instrument (compare with Fig.|2]i. Errors are given as one sigma. 
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To calculate the properties of the lens systems, i.e., masses 
and SFRs, we use available observations in the HST/ACS optical 
filters (f435w, f475w, f555w, f625w, f775w, f8501p), as well as 
HST/NICMOS IR filters (fl lOw, fl60w) from Richard et al. (in 
preparation). 

Using the cluster model to calculate magnifications, we can 
now investigate the observability of such systems. In order to 
be observable, not only must the time delay be smaller than the 
survey length (here assumed to be 5 years) but also the lumi- 
nosity of both images must he above the observation threshold. 
Here we have chosen a magnitude limit of 25.25 mag (Vega) in 
F140W, corresponding to what could be achieved, e.g., with 1/3 
orbit with the HST/WFC3. In Fig. |2] the predicted light curves 
for SNe of Type la and IIP are shown. As can be seen, nearly 
all SNe la exploding in one of the known multiply lensed back- 
ground galaxies would be bright enough to be observable, while 
about half of the SNe IIP, mostly at lower redshift and with 
high magnification, would lie above the luminosity threshold. 
Columns 8-14 in Table |3] indicate which SN type would be ob- 



servable for each image combination producing appropriate time 
delays. 

4. Cluster model constraints 

When modeling the potential of a cluster as described in Sect.|2] 
information on the position of the multiple images is used to 
optimize the free parameters. In general, the addition of a con- 
stant surface density to a lens potential leaves the resulting im- 
age positions of a system unchanged. This is known as the mass 
sheet degeneracy. In the case of A 1689, due to its large number 
of known multiply imaged systems at different redshifts, this de- 
generacy is broken. However, for galaxy clusters with only very 
few observed lens systems, this degeneracy can be problematic. 
The detection of a SN la behind a cluster offers a direct way of 
constraining the len s potential. Due to the standa rd candle na- 
ture of SNe la, see iGoobar & Le ibundgull (1201 ll) for a recent 
review, the information on the absolute magnitude can be used 
to measure the magnification at the position of the SN and com- 
pare it to mass model predictions. One would then be able to 
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Fig. 2: Simulated light curves of SNe la and IIP for pairs of images of strongly lensed galaxies behind A1689. Only pairs with 
time-delay < 5 years are shown. Additional information on different image pairs labeled by, e.g., 1.2 and 1.1 in the top right corner 
of each sub panel, can be found in Table[3] The dashed line indicates an expected magnitude limit 25.25 (Vega) in F140W. 
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rule out realizations producing deviant magnifications. The pre- 
cision of the magnification measurements are ultimately limited 
by the intrinsic scatter in the brightness of SNe la after correc- 
tions for lightcurve shape and color, about 0. 1 mag in the rest- 
frame optical wavelength region (Conlev et al. 201 1). We do not 
anticipate any additional sources of eiTor due to the cluster in 
the line of sight. E.g., observing at near-IR wavelengths ensures 
that corrections from dust in the (low-z) cluster are small, espe- 
cially since galaxy clusters are relatively dust-free environments, 
see Dawson et al. (2009) and references therein. Thus, a SN la 
exploding in any of the background galaxies behind the clus- 
ter could be used to put constraints on the lensing potential. As 
discussed in Paper II, we expect to detect ~ 20-30 SNe la (de- 
pending on the underlying rates estimates) to be detectable, e.g., 
in a 5 year monthly survey at VLT. In order to assess the power 



of this method, we investigate the strength of the correlation be- 
tween the optimized model parameters and the magnification for 
different positions behind the cluster. 

Figure |3] shows the absolute value of the correlation coeffi- 
cients for the predicted magnification and the free input param- 
eters as a function of position for an image at redshift z - 2. 
The parameters shown are those describing the two dark mat- 
ter clumps (Clump 1 and Clump 2) and the cluster galaxy scal- 
ing relations (L*rcut and L*cr) in the model of A1689. A value 
close to 1 implies a very strong correlation while a correlation 
coefficient close to implies no correlation. As can be seen in 
Fig.|3] the free parameters differ in the strength of their coiTela- 
tion with the predicted magnification both with each other and 
the position behind the cluster. The dependence of the correla- 
tion strength on redshift is very weak. In general, of all 33 free 
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Fig. 3: Correlations between the predicted magnification as a function of image position for a source at redshift z — 1 and the input 
parameters which are optimized in the model for A1689 (compare Table [U in a field of view + 100 arcseconds around the cluster 
center Correlations are given as the absolute values of the correlation coefficients. A high value (close to 1; dark red) indicates that 
that the magnification for an image at the given position will display large variations when the input parameter is varied. Being able 
to measure the absolute magification at this image position, using e.g. SNe la observations, will thus make it possible to improve the 
constraints on the cluster parameter If the correlation is low (close to 0; light yellow), we do not expect magnification information 
to significantly improve the current constraint on the cluster parameter at hand. The strongest correlations and best possibility of 
improving the model once a SN la is observed, are the ellipticity, e, core radius, rcore, and velocity dispersion, cr, of the main lensing 
potential. The last panel shows an ACS image of the cluster overlaid with the critical lines for z - 2. The FOV in all panels is 
identical with Fig.[T| 
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Fig. 4: Possible constraints on the model parameters describing the velocity dispersion (left) and core radius (right) of the main dark 
matter clump of A 1689 using information on the magnification from a SN la observed in the multiply lensed background galaxy 
image 18.1 at z = 1.82. The shaded bars represent the parameter distribution from all realizations while the unfilled bars give the 
parameter distribution considering only reaUzations which reproduce the observed magnification within + 0.1 mag. 
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optimized parameters, the ones which show the strongest cor- 
relations and therefore the possibility of improving the model 
once a SN la is observed, are the ellipticity, e, core radius, rcoie, 
and velocity dispersion, cr, of the main lensing potential (Clump 
1). Other parameters, like the corresponding parameters for the 
cluster galaxies and the position of the main lensing potentials 
show a weaker correlation and can probably not be further con- 
strained using SN la observations in the case of A1689. 

Figure |4] shows an example of how the model parameters 
could be constrained if the magnification of an image is known 
to a precision of 0.1 mag. In this example we asssume a SN la 
observed in the multiply-lensed image 18.1 with spectroscopic 
redshift z - 1.82 under the most favorable conditions. Since this 
image does not have a counter image with a time delay below 
five years, it is not included in Table |3] At the position of this 
image, there are strong correlations between the model param- 
eters and the predicted magnification of this image. The corre- 
lation coefficient for the velocity dispersion and core radius of 
the main clump (Clump 1) are 0.84 and 0.87, respectively. As 
can be seen in Figure |4l the possible constraints on cr for the 
main lensing potential improve from ~ 1% to ~ 0.7%, while the 
coiTesponding numbers for rcoie are ~ 4.5% and ~ 2.6%. Thus, 
adding only one constraint from the observed magnification of a 
SN la, the constraints on the model parameters for this already 
very well constrained cluster model of A1689 (1 14 multiple im- 
ages and 24 spectroscopic redshifts), improve by almost a factor 
of 2. It should however be noted, that the possible constraints 
from this image should be seen as a best case scenario, as the 
correlations for other systems are typically weaker, with median 
absolute correlation coefficient values in a FOV extending +100 
arcseconds around the center of the cluster at z = 2 of 0.67 and 
0.7 1 for the velocity dispersion of the main clump (Clump 1) and 
its core radius, respectively. Since we are expecting ~ 5 SNe la 
exploding in all background galaxies per year, we would be able 
to obtain a combined constraint from such observations of the 
magnification at several positions behind the cluster. 

This technique might prove to be even more interesting for 
testing the mass modeling of other less constrained galaxy clus- 



ters. To evaluate the power of the method, we investigate clus- 
ter A2204 with only one known multiply lensed image system, 
making it much less constrained than A 1689, which we have 
been focusing on so far. Again, we investigate the coiTelations of 
the optimized free parameters with the resulting magnifications. 
The correlation strengths for a grid of image positions ±100 
arcseconds around the cluster center at redshift z = 1 is shown 
in Fig. |5] Similarly to the case of A1689, the strongest corre- 
lations are found for cr and rcoie of the main potential with the 
median absolute correlation coefficient values in the FOV shown 
in Fig. |5]at redshift z = 1 being 0.96 and 0.87, respectively. It 
is notable that these correlations are rather weak in the central 
regions of the cluster potential, which is constrained through the 
location of the known multiple images. 

In Figure |6] we show the constraints on these parameters 
from the magnification of a hypothetical SN la with z = 1 close 
to the position of the known lensed image 1.1 measured with 
an uncertainty of + 0.1 mag. The values of the correlation co- 
efficients for cr and rcore of the main clump (Clump 1) are 0.90 
and 0.79, respectively. In this case, the constraints on the velocity 
dispersion for the main lensing potential improve from ~ 17% to 
~ 3.5%, while the corresponding constraint for the core radius 
goes from essentially unconstrained (with a forced upper limit 
from observations) to ~ 25%. Thus, observations of SNe la be- 
hind galaxy clusters promise to be a powerful tool in constrain- 
ing their overall lens potential and thereby the properties of the 
dark matter component. 

5. Measuring the Hubble constant 

It has earlier been suggested to use time delays from multiply- 
lense d sources to measure the Hubble constant (e.g., iRefsdall 
1964'). For a strong lens system, the arrival time of an image 
at angular position 9 relative to the unlensed case for a corre- 
sponding source position /3 is given by 



At : 



1 + zl DlDs 
c Dls 



■m 



(4) 
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Fig. 5: Correlations between the predicted magnification of an image as a function of image position for a source at redshift z - I 
and the different input parameters which are optimized in the model for A2204 (compare Table |2] and Fig. [3]) in a field of view + 
100 arcseconds around the cluster center Correlations are given as the absolute values of the correlation coefficients. The last panel 
shows an ACS image of the cluster overlaid with the critical lines for z = 1 . 
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Fig. 6: Possible constraints on the model parameters describing the velocity dispersion (left) and core radius (right) of the main dark 
matter clump for A2204 using information on the magnification from a SN la observed in the image. The shaded bars represent 
the parameter distribution from all realizations while the unfilled bars give the parameter distribution considering only reaUzations 
which reproduce the observed magnification within + 0.1 mag. 
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where zl is the redshift of the lens, Dl, Ds, and Dls are, respec- 
tively, the angular diameter distances to the lens, to the source, 
and from the lens to the source, and t//{6) is the lens potential. 

Since the ratio of the angular diameter distances scales in- 
versely with the Hubble constant, D = oc 7/^'. By model- 
ing the lens potential, i(/{6), and the source position, j8, one can 



use time delay measurements from strongly lensed systems to 
estimate the Hubble constant. 

There is also a dependence on D from other cosmological pa- 
rameters, such as Q.M and Qa, although this de pendence is much 
weaker than that on the Hubble constant (e.g.. | Bolton & BurlesI 
l2003HCoe & Moustakasll2009l:ISuvu et al]l20m 
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Fig. 7: Predicted time delays in days and percentual one sigma 
errors for all the image pairs of multiply imaged background 
galaxies with spectroscopic redshift and time delay smaller than 
five years behind A1689. The color of the markers indicates the 
redshift of the systems. 
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Fig. 8: Constraints on the time delay between images 19.3 and 
19.4 at z = 2.60 behind A1689 using the magnification of both 
images. In the shaded histogram the time delay distribution from 
the complete Monte Carlo Markov Chain is shown. The unfilled 
histogram shows the distribution after the cuts due to magnifica- 
tion constraints. Errors are given as one sigma. 
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Due to the transient nature of SNe, the observed time de- 
lay from SN light curves can be determined with high precision, 
typically on the order of less than a few days. To improve the 
constraints on //(> it is thus important to be able to narrow down 
the errors of the lens potential. Here we investigate the con- 
straints on //() possible when observing a SN in one of the known 
multiply-imaged background galaxies of A 1689 with predicted 
time delays below 5 years and spectroscopic redshifts. Summing 
up our SN rate estimates (given in Table |3]l, we expect an aver- 
age total of 3 SNe exploding in those galaxies during a survey 
time of 5 years. The effective number of strongly lensed SNe 
which are observable, will depend on the survey properties. For 
a five year monthly survey using the HST/WFC3, we would ex- 
pect to detect approximately 1 strongly lensed SN in these galax- 
ies. However, one should keep in mind that there are 10 more 
known strong lensing systems behind A 1689, consisting of 30 
images, which have not been included in this study since they 
have photometric redshifts only. Furthermore, as observations of 
this massive cluster c ontinue, additional multip ly lensed systems 
might be discovered dCoe & MoustakaJ2009l) . 

From the different realizations of the cluster modeling, we 
get a prediction with error estimates for the time delays between 
image pairs in the known multiply lensed background galaxies. 
These errors depend on the positions of the image pairs. As can 
be seen in Fig.|7] for more than half of the image pairs with pre- 
dicted time delays below five years and spectroscopic redshift, 
the corresponding errors are below 20%. For several systems the 
time delay precision from the model lies around 10% and below. 
Eventually combining constraints from several strongly lensed 
SNe, would allow the value of the Hubble constant to be deter- 
mined with high precision. Thus, nominally, to match the current 
~ 3% accuracy of the local measurements of //q (Riessetal. 
1201 Ih . about ten massive clusters should be monitored for five 
years. Although it may seem like a rather large time investment 
in telescope time to match a result already at hand, we empha- 
size the importance of a non-local test of the universal Hubble 
scale, a key test in cosmology. 

In the case of a lensed SN la, we can use the information 
on the absolute magnification to rule out realizations that predict 



magnifications outside the range allowed by observations. With 
this extra information, we might be able to further constrain the 
time delay error from the lens model. Therefore, we check for 
correlations between the predicted magnification of the images 
in a system and the time delay between the images. While many 
systems show no or only weak correlations between magnifica- 
tions and time delay, there are a few more promising systems. 

Figure [8] shows an example for the cuts in time delay er- 
ror which can be made using magnification constraints from a 
lensed SN la in background galaxy 19. In this case, the errors 
on the time delay between images 19.3 and 19.4 can be reduced 
from ai^jl^ to %. 

As discussed in Sect. |4] the additional model constraints 
from a measurement of the magnification by using the standard 
candle nature of SNe la are only modest for a well constrained 
cluster as A1689. However, such a constraint would prove much 
more valuable in the case of a time delay measured in a less 
constrained cluster as it could break t he degeneracy between th e 
mass profile and the Hubble constant (lOguri & Kawandl2003h . 

6. Discussion and Conclusions 

A SN la exploding in any of the background galaxies of a clus- 
ter will give an absolute measure of the magnification at that 
point and can therefore potentially be used to constrain the mass 
modeling of the cluster For cluster A 1689, which we have been 
focusing on in this study, the cluster potential is already very 
well modelled. However, such an additional constraint from one 
or more lensed SNe could help to improve the cluster model fur- 
ther by reducing parameter uncertainties. For less constrained 
galaxy clusters, this technique proves to be even more interest- 
ing. Using the magnification information from a single SN la 
behind cluster A2204 can be a powerful tool to narrow down the 
model parameters of the dark matter halo component. 

Since SNe are point-like sources for a limited period, mi- 
crolensing from stars in a cluster galaxy close to the line of 
sight might in principle have an important effect on the mag- 
nification and thus affect the possible constraints on the large 
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scale mass model from S Ne la. However, it has been concluded 
(lOguri & Kawano"2003) that this effect should only cause small 
deviations. Microlensing events caused by massive compact halo 
objects (MACHOs) in the intracluster medium might unambigu- 
ously be identified in SN light curves. Thus any modulation of 
the light curve shape, or lack thereof, can put limits on the mass 
fraction of the cluster mass in the form of 10"^ < M/Mq < 10""* 
MACHOs (Kolatt & Bartelmann 1998). 

In the case of a SN exploding in one of the known multi- 
ply lensed background galaxies, there is a good chance of ob- 
serving multiple images of the SN and determine the observed 
time delay with high precision. The intrinsically brighter SNe, 
e.g. SNe la and SNe lln, should in principle be detectable in 
essentially all of the systems. However, for SNe IIP, which are 
expected to be the most common SN type, the magnification is 
only strong enough for them to become observable in about half 
of the cases. 

Due to the well constrained mass model for galaxy cluster 
A 1689, SN time delays can be used as an independent measure 
on the Hubble parameter at high redshift. In case of a SN la, the 
magnification might also be used as a constraint to further im- 
prove the error in the time delay predicted by the lensing model. 
Another advantage with using a cluster with a well constrained 
mass model is the possibility of predicting subsequent images, if 
a SN behind the cluster is observed. 

It is quite exciting that CLASH, a clu ster monitoring multi - 
cycle program is currently been pursued (iPostman et al.ll20TTl) . 
However, the strategy in the CLASH program is tuned for find- 
ing and accurately studying SNe la in the parallel fields, i.e., in 
the low magnification region. Much of the observations of the 
cluster cores are done in optical and UV filters, where high-z 
SNe would not be detectable. Furthermore, the cluster fields are 
only monitored for 2-4 months making the search for multiple 
images practically impossible, especially considering the time 
dilation (1 + z) =2-4.5 in the lensed SN lightcurves. 

Corroborating the value of Hq found in the local universe 
could provide crucial support to our cosmological picture. On 
the other hand, a discrepancy would falsify the accepted scenario 
and shed new light into the dark matter and dark energy puzzles. 
Since strongly lensed SNe involve both angular diameter and lu- 
minosity distances {(Ja and c/l), such systems provide an unique 
chance to probe the distance reciprocity relation, c/l - i^+z)^dA- 
Any violation of this relation would indicate either the existence 
of unaccounted astrophysical systematic effect affecting one of 
the distance measure but not the other, or that the Universe is not 
described by the standard cosmological model. Specifically, a 
violation could indicate that the Universe is described by a non- 
metric theory of gravity or a theory in which light does not travel 
on unique null geodesies. Curren t tests of di^/dA are restricted to 
z < 1 (e.g.. lLampeitl et al]|2010t) and have > 20% uncertainties. 
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